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Abstract. Nonlinear self-modulation of the dust acoustic waves is studied, in the 
presence of non-thermal (non-Maxwellian) ion and electron populations. By em- 
ploying a multiple scale technique, a nonlinear Schrodinger-type equation (NLSE) 
is derived for the wave amplitude. The influence of non-thermality, in addition to 
obliqueness (between the propagation and modulation directions), on the condi- 
tions for modulational instability to occur is discussed. Different types of localized 
solutions (envelope excitations) which may possibly occur are discussed, and the 
dependence of their characteristics on physical parameters is traced. The ion de- 
viation from a Maxwellian distribution comes out to be more important than the 
electron analogous deviation alone. Both yield a de-stabilizing effect on (the am- 
plitude of) DAWs propagating in a dusty plasma with negative dust grains, and 
thus favor the formation of bright- (rather than dark-) type envelope structures 
(solitons) in the plasma. A similar tendency towards amplitude de-stabilization is 
found for the case of positively charged dust presence in the plasma. 



1. Introduction 

Dust contaminated plasmas have recently received considerable interest due to their 
wide occurrence in real charged particle systems, both in space and laboratory 
plasmas, and to the novel physics involved in their description [1, 2]. An issue of 
particular interest is the existence of special acoustic-like oscillatory modes, e.g. 
the dust-acoustic wave (DAW) and dust-ion-acoustic wave (DIAW), which were 
theoretically predicted about a decade ago [3, 4] and later experimentally confirmed 
[5, 6]. The DAW, which we consider herein, is a fundamental electrostatic mode in 
an unmagnetized dusty plasma (DP), which has no analogue in an ordinary e-i 
(electron - ion) plasma. It represents electrostatic oscillations of mesoscopic size, 
massive, charged dust grains against a background of electrons and ions which, 
given the low frequency of interest, are practically in a thermalized (i.e. Boltzmann) 
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equilibrium state. The phase speed of the DAW is much smaller than the electron 
and ion thermal speeds, and the DAW frequency is below the dust plasma frequency. 

The amplitude (self-) modulation of the DAse waves has been recently considered 
[7 - 9] , by means of the reductive perturbation formalism [10] . A study of the mod- 
ulational stability profile of DA waves has shown that long/short wavelength DA 
waves are stable/unstable against external perturbations. The respective parameter 
regions are associated with the occurrence of localized envelope excitations of the 
dark/bright type (i.e. voids/pulses). Obliqueness in perturbations are found shown 
to modify this dynamic behaviour [7 - 9] . 

Allowing for a departure from Boltzmann's distribution for the electrostatic back- 
ground {non-thermality) has been shown to bear a considerable effect on electro- 
static plasma modes. Inspired by earlier works on the ion-acoustic (IA) solitary 
waves [11], recent studies have shown that the presence of a non-thermal ion ve- 
locity distribution may modify the nonlinear behaviour of the DA waves, affecting 
both the form and the conditions for the occurrence of the DA solitons [12, 13]. 
Also, the self-modulation of IA waves was recently shown to be affected by the 
electron non-thermality [14]. However, no study has been carried out of the effect 
of a non-thermal ion and/or electron background population on the modulation 
properties of the DA waves. This paper aims at filling this gap. 



2. The model 

Let us consider the propagation of the dust-acoustic waves in an unmagnetized 
dusty plasma. The mass and charge of dust grains (both assumed constant, for 
simplicity) will be denoted by m& and qd = s Z^e, where s = sgnq^ = qd/\<ld\ 
denotes the sign of the dust charge (= ±1). Similarly, we have: m, and % = +Z,fi 
for ions, and m e and q e = —e for electrons. The typical DAW frequency will be 
assumed much lower than the (inverse of the) characteristic dust grain charging 
time scale, so charge variation effects (known to lead to collisionless wave damping) 
will be omitted here. 

The basis of our study includes the moment - Poisson system of equations for the 
dust particles and non-Boltzmann distributed electrons and ions. The dust number 
density rid is governed by the continuity equation 

^+W-(n d u d ) = 0, (2.1) 
and the dust mean velocity obeys 

^+u d .Vu d = -^V$ (2.2) 
at nid 

where $ is the electric potential. The dust pressure dynamics (i.e. the dust temper- 
ature effect) is omitted within the cold dust fluid description. The system is closed 
by Poisson's equation 

V 2 $ = -4tt e (m Z l -n e + s n d Z d ) . (2.3) 

Overall neutrality, viz. 

n ifi Zi - n efl + s n d .o Z d = 



is assumed at equilibrium. 
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2.1. Non-thermal background distribution(s) 

Following the model of Cairns et al. [11], the non-thermal velocity distribution for 
e~ (electrons) and i + (ions) is taken to be 

, / s »V,o 1 + a s'V i /vf h , 

f s >{v;a s ,) = — — exp(-t> /2v th ,) , (2.4) 

yj^ks' 1 + 3as ' 

where ?v ; o, TV and u^g' = (fesTV/TTv) 1 / 2 denote the equilibrium density, tem- 
perature and thermal speed of the species s' € {1,2} = {i, e}, respectively. The real 
parameter a s > expresses the deviation from the Maxwellian state (which is recovered 
for <v = 0). Integrating over velocities, Eq. (2.4) leads to 

n.> = n a /, (l + a s >b s >& + b s ,<S>' 2 ) exp(-<7 s >$') (2.5) 

[11], where <J\/2 — o~i/ e = +1/ — 1 and b s > = Aa s >/(1 + 3a s >); the normalized 
potential variables are $' s , = Z s ie$>/kBT s >, where Z 1 / 2 = Zi/1 is the ion/electron 
charge state. 

2.2. Reduced model equations - weakly nonlinear oscillations 

The above system of evolution equations for the dust fluid can be cast into the 
reduced (dimensionless) form: 

On 

— +V(nu) = 0, (2.6) 

and 

-^T+u- Vu = -sVcj) (2.7) 

where the particle density rid, mean fluid velocity Ud, and electric potential $ are 
scaled as: n — rid/ndfi, u = Ud/cd, and <fr — \qd\^ / {ksT e f f), where nd.o is the 
equilibrium dust density; the effective temperature T" e // is related to the char- 
acteristic dust speed Cd = (ksTeff/md) 1 ^ 2 — uj p .d\D,eff , defined below (fee is 
Boltzmann's constant). Time and space variables are scaled over the characteris- 
tic dust period (inverse dust plasma frequency) oj~ d = ('lnndfiqd/ m d)~ 1 ^ 2 , and 
the effective Debye length, defined as \D, e ff = [(1 — &e)A^ e + (1 — b^X^ 2 ^ 1 ^ 2 
[where Xn.e/i — (kBT e / i / / i'Kn e / ifi q 2 / i ) 1 / 2 is the Debye length for species e/i; the 
non-thermality parameters b e j i are defined below]. The (unperturbed) dust Debye 

length \d4 = (^rfe + ^rfi) 1 ^ 2 ^ a ^ so defined. 

Near equilibrium (<j> <C 1), Poisson's Eq. (2.3) becomes: 

V 2 </> w <j)-acj) 2 + a'(t> 3 -s(n-l). (2.8) 

Note that the right-hand side in Eq. (2.8) cancels at equilibrium. Here, 

Zi 1 \ X D,eff T eff 



and 



X D,i T i X D, e Te J 2Z d 



2 1 , iu 1 \ 2 „ tf T 2 



{l + 3bi)Zf , 1 + 36, 



+ 



\2 rp2 1 \2 rp2 

^D.i 1 ! A D,e 1 e 



Ml 



For T e 3> Ti, one may retain the approximate expressions: a w — s(l — /i)/[2(l — 6 2 ) 2 ] 
and a' w (1 + 36;)(1 - /i) 2 /[6(l - & 2 ) 3 ]; also, A Ae// w A Ai (l - b t )-^ 2 in this 
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case (notice that the dependence on the electron parameters disappears in this 
approximation). We have defined the dust parameter \i = n e ,o/(Zin>i,o) = 1 + 
sZdridfil '(ZiHifi); see that fj, is lower (higher) than unity for negative (positive) 
dust charge, i.e. for s = — 1 (+1). 

The dust parameter ji, together with the temperature ratio U = Ti/T e , are the 
(dimcnsionlcss) physical parameters essentially affecting our problem. For instance, 
see that Xn,d (defined above) may be expressed as: Xda = Xo, e {(ti/j,/Zi)/[l + 



The system of the evolution equations (2.6) - (2.8) will be the basis of the analysis 
that follows. 



3. Perturbative analysis 

Following the reductive perturbation technique [10], we define the state vector 
S = {n, u, <f>}, whose evolution is governed by Eqs. (2.6) - (2.8), and then ex- 
pand all of its elements in the vicinity of the equilibrium state S^ ' = (1, 0, 0) T , viz. 
S = S(°) +eSW +e 2 S( 2 ' + ..., where e <C 1 is a smallness parameter. We assume that 



Sj (n) = J2Z-oo S$(X, T)e a( ^ r -^ (for j = 1,2,...; we impose S™, = 5] 



coordinates X = e(x — Xt) and T — e 2 t, where A (G 3?) (bearing dimensions of ve- 
locity) will be determined later. Note that the choice of direction of the propagation 
remains arbitrary, yet modulation is allowed to take place in an oblique direction, 
characterized by a pitch angle 9. Having assumed the modulation direction to define 
the x— axis, the wave vector k is taken to be k = (k x , k y ) — (k cosO, k sin#). 

According to the above considerations, the derivative operators in the above 
equations are treated as 



(W^)]} 1/2 - 




d d d 2 d 



V 



V + ex — 



and 



V 2 -» V 2 + 2e 



dxd( 




2 ' 



i.e. explicitly 





and 




for any 



A\ n> of the components of S\ 
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3.1. First harmonic amplitudes 

Iterating the standard perturbation scheme [10], we obtain the first harmonic am- 
plitudes (to order ~ e 1 ) 



s(l + k 2 )<P, 



,(!) - 



<i.x = ^costtnl 



(1) w • a (!) 
u\ y = — sint/nj 



(3.1) 



which may be expressed in terms of a sole state variable, e.g. the potential correction 
<f>^ = ip. The well-known DA wave (reduced) dispersion relation lu 2 = k 2 / (k 2 + 1) 
is also obtained, as a compatibility condition. 

The amplitudes of the 2nd and 0th (constant) harmonic corrections are obtained 
in order ~ e 2 (the lengthy expressions are omitted for brevity). Furthermore, the 
condition for suppression of the secular terms leads to the compatibility condition: 
A = du/dk x — uj'(k)cos9 — kcos9/[uj(l + fc 2 ) 2 ]; therefore, A bears the physical 
meaning of a group velocity in the direction of amplitude modulation (<~ x) . 

3.2. The envelope evolution equation 

The potential correction tp obeys an explicit compatibility condition in the form of 
the nonlinear Schrddinger-type equation (NLSE) 



dip d 2 ^ 
8T + OX* 



Q\iP\ 2 iP = 0. 



(3.2) 



The dispersion coefficient P is related to the curvature of the dispersion curve as 
P = d 2 oj/2dk 2 x = [u"{k) cos 2 + u'(k) sin 2 6/k] /2; the exact form of P reads 



m = 2^(i) [i-(i+^- 2 )«»« 2 »- 



(3.3) 



The nonlinearity coefficient Q, due to carrier wave self-interaction, is given by a 
complex function of k, a and a' . Distinguishing different contributions, Q can be 
split into three distinct parts, viz. 



where 



Qo = + 



1 



1 



Q = Qa + Qi + Qi ■ 



l 



(3.4) 



2co (1 + fc 2 ) 2 1 - A 2 
x\k 2 [3 + 6fc 2 + 4fc 4 + fc 6 + 2a(s(2fc 2 + 3)+ 2a A 2 )] 

+ (2 + 4fc 2 + 3fc 4 + fc 6 + 2sa(3) cos 29 
+ 2(l + fc 2 ) 4 ^ 2 cos 2 6 



Qi 
Q 2 



+ fc(l + fc 2 ) 

3 a' k 2 



k 2 +uj 2 (l + k 2 ) 



- (1 + k 2 + 2sa(3) cos6 



2uj (1 + fc 2 ) 2 ' 
1 1 1 



(3.5) 
(3.6) 



12 w fc 2 (l + fc 2 ) 2 
x 1 2 k 2 [5 s a (1 + fc 2 ) 2 + 2a 2 + 3(1 + k 2 f (1 + 3fc 2 )] 
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+ (1 + fc 2 ) 3 u 2 (3 + 9fc 2 + 6fc 4 + 2sa) j . (3.7) 

One may readily check, yet after a tedious calculation, that expressions (3.5) and 
(3.7) reduce to (53) and (54) in Ref. 7 (for a Maxwellian distribution, i.e. setting 
o,i = ci2 = in all formulae above); notice however, that the term Qi was absent 
therein. Note that only Q (which is related to self-interactions due to the zeroth 
harmonic) depends on the angle 9 (and is, in fact, an even, it— periodic function of 
9). 

4. Stability profile envelope excitations 

The evolution of a wave whose amplitude obeys Eq. (3.2) essentially depends on 
the sign of the product PQ (see, e.g., in Ref. 16 for details; also see in Refs. 9 
and 17) which may be numerically investigated in terms of the physical parameters 
involved. 

For PQ > 0, the DA wave is modulationally unstable and may either collapse, 
when subject to external perturbations, or evolve into a series of bright-type local- 
ized envelope wavepackets, which represent a pwZse-shaped solution of the NLSE 
(3.2). This type of solution is depicted in Fig. 2. 

For PQ < 0, on the other hand, the DA wave is stable and may propagate as 
a dark/ grey '-type envelope wavepackct, i.e. a propagating localized envelope hole 
(a void) amidst a uniform wave energy region. Notice the qualitative difference 
between the dark and grey solutions (depicted in Fig. 3a and b, respectively): the 
potential vanishes in the former, while it remains finite in the latter. 

The exact analytical expressions for different types of envelope excitations (de- 
picted in Figs. 2 and 3) are summarized in Refs. 9 and 17; more details can be 
found in Ref. 18. We note that, in either case (i.e. for positive or negative prod- 
uct PQ), the localized excitation width L and maximum amplitude ipo satisfy 
Lipo — (2P/Q) 1 / 2 = constant. The dynamics, therefore, essentially depends on 
the quantity n = P/Q, whose sign (magnitude) will determine the wave's stability 
profile and the type (characteristics, respectively) of the localized envelope excita- 
tions which may occur. 

5. Numerical analysis 

According to the above analysis, both the DAW stability profile and the type of DA 
envelope excitations possibly occurring in the plasma are determined by the sign 
of the product of the NLSE coefficients P and Q, which is essentially a function 
of the wavenumber k (normalized by kjj = A^, 1 ^ [19]), the dust parameter fi, the 
temperature ratio Ti/T e = ti and the non-thermality parameters dj ;e , in addition 
to the modulation angle 9. The exact expressions obtained for P and Q may now 
be used for the numerical investigation of the wave's modulational stability profile. 

In the figures provided here (see Figs. 4 - 9), the black/white regions represent a 
negative/positive PQ sign, implying DAW being modulationally stable/ unstable 
(and possibly propagating in the form of dark/bright type excitations, respectively). 

Throughout this Section, we have used the set of representative values: ti = 
Ti/T e = 0.1, Z, = 1 and 5 = Z d n dfi / (Zin ifi ) = 0.25, implying fj, = 0.75 for 
negatively charged dust (and /i = 1.25 for positively dust charge; cf. the definitions 
above). The (normalized) wavenumber fc, modulation angle (0 < 9 < n) and the 
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non-thermality parameters (0 ^ a e ,i < 1) will be allowed to vary according to the 
focus of our investigation. The negative dust charge case (s = —1) is implied, unless 
otherwise mentioned. 

5.1. Parallel modulation 

In the special case of the parallel modulation (6* = 0), the analytical behaviour of P 
and Q was briefly studied in Ref. 15. See that P = -3uj 5 /{2k i ) < for 8 = so that 
one only has to study the sign of Q. For k <C 1 one then obtains the approximate 
expression [15]: P ~ — 3fc/2 and Q w (3 + 2sa) 2 /(\2k) > 0, which prescribes 
stability (and dark/grey type envelope excitations - voids) for long wavelengths. 
For shorter wavelengths, i.e. for k > k cr (where k cr is some critical wavenumber, 
to be determined numerically), Q changes sign (becomes negative), and the wave 
may become unstable. 

In addition to the above results, an increase in the non-thermal ion population 
(i.e. in the value of at) is seen to favor instability at lower values of k: see that 
the black region shrinks - for 8 = 0- when passing from Fig. 4a to Fig. 4b. Large 
wavelengths always remain stable and may only propagate as dark-type excitations. 
On the other hand, the effect of an increase in the non-thermal electron population 
(i.e. in the value of a e ), while keeping ions Maxwellian, has a similiar, yet not so 
dramatic effect; cf. Fig. 4c, for 8 = 0. We note that the effect of non-thcrmality 
(of ions, in particular) on the parallel modulation of a dusty plasma with positively 
charged dust is qualitatively similar; compare e.g. Fig. 4b to Fig. 7b (for 8 = 0). 

5.2. Obliqueness effects 

Let us now consider an arbitrary value of 8. First, we note that, for small wavenum- 
ber k, one obtains: P w sm 2 8/(2k) > and Q w -(3 + 2sa) 2 /(l2k) < 0, so that 
stability (and the existence of dark/grey type excitations) is, again, ensured for 
long wavelengths A 3> Xd- 

We proceed by a numerical investigation of the sign of the coefficient product 
PQ in the (k, 8) plane; see Fig. 4. The profile thus obtained is qualitatively rem- 
iniscent of the results in Ref. 9: the critical wavenumber, say k cr , beyond which 
PQ changes sign, decreases as the modulation obliqueness angle 8 increases from 
zero (parralel modulation) to - roughly - 0.5 rad (i.e. approximately 30), implying a 
destabilization of the DAW amplitude for lower wave number (higher wavelength) 
values. Now, above 0.5 rad (i.e. 60), approximately, and up to n/2 rad (i.e. 90, 
transverse modulation) the wave remains stable. 

A question naturally arises at this stage: how does the forementioned stability 
profile depend on the non thermal character of the background? The ion and elec- 
tron non-thermality effects are separately treated in the following paragraphs. 

5.3. Ion non-thermality effects 

Let us vary the value of the ion non-thermality parameter Oj. Passing to a finite 
value of di = 0.2 results in a narrower stability region for small angles (i.e. k cri 
at which PQ(k cr ) = 0, decreases with o», for small 8) - see Fig. 4b - while the 
influence on the behaviour for larger angles is less important. Therefore, a small non 
Maxwellian ion population (for tij less than 0.2, roughly) seems to slightly enhance 
modulational instability (i.e. destroy stability) at lower wavenumber values, at least 
for a small to moderate modulation obliqueness. Slightly higher values of dj rather 
favor instability even further. Finally, the picture is reversed for very high cij (above, 
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say, 0.5), where the wave is slightly stabilized by increasing a,i (at least for moderate 
values of 9). These effects may be depicted by keeping the value of 9 fixed, while 
varying the value of af, see Figs. 5. The influence of non-thermality on strongly 
oblique DAW modulation is tacitly negative: stability is destroyed by increasing at 
(see e.g. Fig. 5c, where 9 — 7r/4). The transverse DAW modulation case (i.e. for 
9 = it/2) remains globally stable, so the corresponding (black) diagram was omitted 
in Fig. 5. 

5.4. Electron non-thermality effects 

In analogy to the previous paragraph, we may now vary the value of the electron 
non-thermality parameter a e . The effect is qualitatively similar, yet far less dra- 
matic; see Fig. 6. This is physically expected, since the influence of the background 
ions on the inertial dust grains is more important than that of the lighter electrons. 

5.5. Positive dust 

The analogous investigation in the positive dust case (s = +1; /1 = 1.25 everywhere 
below) reveals a qualitatively different picture. Again, the addition of positive dust 
is seen to result in a much narrower stability region; compare Figs. 4 and 7. We 
see that positive dust slightly favors instability, with respect to negative dust; cf. 
the discussion in Ref. 9. 

As far as the influence of ion non-thermality on stability is concerned, the qual- 
itative effect seems to be the opposite as in the preceding paragraphs. For a given 
value of the angle 9, increasing the ion "non-Maxwellian" parameter a, seems to 
result in a more unstable configuration, i.e. in a narrower black region in Fig. 8, for 
9 = 0; this is true for values of dj up to, say, 0.5, while reaching higher ones bears 
the opposite effect (see the upper half of Fig. 8a). For higher 9, the result is similar, 
yet less dramatic; see Figs 8b, c. Again (as in the positive dust charge case), the 
transverse DAW modulation case (i.e. for 9 = n/2) remains globally stable, so the 
corresponding (completely black) diagram was omitted in Fig. 8. 

Finally, taking into account electron non-thermality, i.e. increasing the electron 
parameter a e from zero to a finite value yields a similar effect, yet much less im- 
portant effect, for practically all values of 9 and k. 

In conclusion, the modulational instability of DAWs propagating in a dusty 
plasma with positively charged dust grains seems to be enhanced by the appearance 
of a non-thermal background component. 



6. Conclusions 

We have investigated the amplitude modulation of dust acoustic waves in the pres- 
ence of a non-thermal (non-Maxwellian) ion and/or electron background, focusing 
on the influence of the latter on the stability profile, and also on the conditions for 
the occurrence of envelope excitations. 

Summarizing our results, we find that the presence of non-Maxwellian electron 
an/or ion distribution(s) is expected to yield a destabilizing effect on the DAW 
propagation in a plasma with negatively charged dust grains (slightly favoring 
dark- rather than bright- type envelope wave packets), in particular for a moder- 
ate obliqueness in amplitude perturbation. This effect is (qualitatively similar but) 
stronger for nonthermal ions rather than electrons. In the presence of a transverse 
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modulation, the wave's stability profile remains intact. The modulational instabil- 
ity of the DAWs propagating in a dusty plasma with positively charged dust grains 
(qd > 0) also seems to be enhanced by the appearance of a nonthermal electron 
and ion components. Again, this effect is stronger for nonthermal ions rather than 
electrons; also, transversely modulated DAWs remain stable. 

Finally, the occurrence of bright (rather than dark) type envelope wave packets 
is rather favored by deviations of the ion (and electron, though less) population(s) 
from the Maxwellian distribution. 

Therefore, adding a non thermal ion and/or electron component to the plasma 
may control (and possibly stabilize) the propagation of dust-acoustic envelope soli- 
tons, by enhancing energy localization via instability of the wave's amplitude (due 
to carrier self-modulation). 
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For rigor, the (reduced) parameter k (i.e. essentially <~ A.o je ///A = k/kD, e ff, 
where A is the wavelength) appearing in the formulae, has had to be redefined 
as k — > k\D, e ff/^D,d, in order to be represented in the horizontal axis in Figs. 
4-9 (where the normalized variable is now <~ An^/A = k/kD,d)- 
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Figure captions 



Figure 1: 

The non-thermal distribution f(v;a), as defined by Eq. (2.4) [scaled by f(v; 0) = 
l/y/2n], vs. the reduced velocity v/v t h, for a = 0, 0.1, 0.2, 0.3 (from top to bottom). 

Figure 2: 

Bright type localized envelope modulated wavepackets (for PQ > 0) for two 
different (arbitrary) sets of parameter values. 

Figure 3: 

(a) Dark and (b) grey type localized envelope modulated wavepackets (for PQ < 
0). See that the amplitude never reaches zero in the latter case. 

Figure 4: 

The sign of the product PQ is depicted vs. (normalized) wavenumbcr k (hori- 
zontal axis) and modulation angle 9 (vertical axis), via a color code: black (white) 
regions denote a negative (positive) sign of PQ, implying stability (instability), and 
suggesting dark (bright) type envelope soliton occurrence. Here, \x = 0.75 (negative 
dust charge) and: (a) a, = a e = (Maxwellian background); (b) = 0.2, a e = 
(ion non-thcrmality) ; (c) a, = 0, a e = 0.2 (electron non-thcrmality) . Remaining 
parameter values as cited in the text. See that the effect in (c) is rather negligible. 

Figure 5: 

The sign of the product PQ is depicted vs. (normalized) wavenumber k (hori- 
zontal axis) and ion non-thermality parameter ai (vertical axis), for a modulation 
angle 9 equal to: (a) (parallel modulation); (b) 0.4; (c) 7r/4; the case 9 = ir/2 - 
transverse modulation - is omitted, since globally stable. Same color code as in Fig. 
4: black (white) regions denote a negative (positive) sign of PQ. Here, \i = 0.75 
(negative dust charge) and a e — (Maxwellian electrons). Remaining parameter 
values as cited in the text. 

Figure 6: 

Similar to the preceding Figure, but for non-thermal electrons and Maxwellian 
ions: the sign of the product PQ is depicted vs. (normalized) wavenumber k (hori- 
zontal axis) and electron non-thermality parameter a e (vertical axis), for a modula- 
tion angle 9 equal to: (a) (parallel modulation); (b) 0.4; (c) 7r/4; the case 9 = ir/2 
- transverse modulation - is omitted, since globally stable. Same color code as in 
Fig. 4: black (white) regions denote a negative (positive) sign of PQ. Here, fi = 0.75 
(negative dust charge) and ai = (Maxwellian ions). Remaining parameter values 
are cited in the text. 



Figure 7: 
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Similar to Figure 4, but for a positive dust charge; here, fi = 1.25. The sign of 
the product PQ is depicted vs. (normalized) wavenumber k (horizontal axis) and 
modulation angle 9 (vertical axis), for: (b) a, = 0.2, a e = (ion non-thermality) ; 
(c) di = 0, a e = 0.2 (electron non-thermality). Remaining parameter values as in 
Fig. 4. See that the effect in (c) is not as dramatic as in (b). 

Figure 8: 

Similar to Figure 5, but for a positive dust charge; here, \i = 1.25 and a e = 
(remaining parameter values as in Fig. 5). The sign of the product PQ is depicted 
vs. (normalized) wavenumber k (horizontal axis) and ion non-thermality parameter 
cii (vertical axis), for a modulation angle 9 equal to: (a) (parallel modulation); (b) 
0.4; (c) 7r/4; the case 9 = tt/2 - transverse modulation - is omitted, since globally 
stable. 



Figure 9: 

Similar to Figure 6, but for a positive dust charge; here, \i = 1.25 and = (re- 
maining parameter values as in Fig. 6). The sign of the product PQ is depicted vs. 
(normalized) wavenumber k (horizontal axis) and electron non-thermality parame- 
ter a e (vertical axis), for a modulation angle 9 equal to: (a) (parallel modulation); 
(b) 0.4; (c) 7r/4; the case 9 = n/2 - transverse modulation - is omitted, since glob- 
ally stable. See that the electron non-thermality effect is less important than its ion 
counterpart (cf. Fig. 6). 
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Figure 8. 
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